Inflationary Super-Hubble Waves and the Size of the Universe 
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The effect of the scalar spectral index on inflationary super- Hubble waves is to amplify/damp large 
wavelengths according to whether the spectrum is red (n s < 1) or blue (n s > f). As a consequence, 
the large-scale temperature correlation function will unavoidably change sign at some angle if our 
spectrum is red, while it will always be positive if it is blue. We show that this inflationary filtering 
property also affects our estimates of the size of the homogeneous patch of the universe through the 
Grishchuk-Zel'dovich effect. Using the recent quadrupole measurement of ESA's Planck mission, we 
find that the homogeneous patch of universe is at least 87 times bigger than our visible universe if 
we accept Planck's best fit value n s — 0.9624. An independent estimation of the size of the universe 
could be used to independently constrain n s , thus narrowing the space of inflationary models. 
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I. INTRODUCTION 

The temperature power spectrum coefficients of the 
Cosmic Microwave Background (CMB) radiation field, 
the Cg&, when plotted as a function of their multipoles 
I, comprise one of the most famous figures of modern 
cosmology. In face of the recent and unprecedented boost 
that this subject has received from Planck's data release 
[1, 2], CMB will certainly retain its distinguished status 
among cosmologists during the next coming years. Not 
so famous, but yet conveying the same information, is the 
plot of the angular two-point correlation function C(9): 



C{9) = (ATxAr 2 ) = J2 



(21 + 1) Ci 



i=i 



4tt 



PeicosO). (1) 



While either the position space C(9) or the harmonic 
Ces are equivalent and sufficient to fully characterize 
the CMB morphology in the Gaussian and statistically 
isotropic regime, the latter are mainly used for extracting 
cosmological information since, in this regime, the Ces de- 
couples from one another, making the harmonic analysis 
much more transparent. The angular correlation function 
in Eq. (1), instead of splitting the information in several 
independent multipoles, goes on the opposite direction 
and compresses information from virtually all multipolar 
scales to give the physics at a single angular separation 0. 
This does not mean that a real space approach to CMB 
is prohibitive, and indeed non-trivial physics can result 
from it [.'!-' ]. However, the result of this compression is a 
featureless and literally monotone function of 6, as shown 
by the continuous line of Fig. (1). Note that, for the 
sake of exposition, we have included a primordial dipolc 
C\ = 3C2 in Eq. (1), with the proportionality constant 
fixed by the Sachs- Wolfe effect £(£ + l)Ce = constant. 

By eye, the most prominent feature of Fig. (1), if not 
the only, is that C(9) has a root at 9 ~ 75°. The imme- 
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Figure 1: CMB two-point correlation function as a function of 
angular separation 9 (solid, thick line) in the ACDM model. 
This curve includes a primordial dipole C\ — 3C2, in accor- 
dance with the Sachs- Wolfe effect. The removal of this dipole 
results in the appearance of two roots, one at ~ 40° and 
other at m 120°; this is shown by the dashed light curve, 
which also includes cosmic variance errors bars. 



diate question which comes to mind is: why is there an 
angular scale above which CMB photons become uncor- 
related? At these angles, the photons we detect are essen- 
tially sourced by plane wave perturbations of the gravi- 
tational potential through the Sachs- Wolfe effect. Since 
these waves have a length comparable to the radius of our 
particle horizon, each of their crests will sub-intend an 
angle of approximately 90° in the CMB sky, while larger 
waves will only have a fraction of their crests inside the 
horizon. Therefore, photons separated by angles larger 
than 90° are, on average, probing the maxima and the 
minima of these waves. This gives us the anti-correlation 
pattern observed, since C(9) = (H — ) < 0. What about 
larger waves which are shifted to the left /right, so that 
only their positive/negative crests lies inside the horizon? 
Photons probing these waves will always be correlated, 
since C(9) = (++) = ( ) > 0. However, since a pri- 
ori there is no upper limit for the wavelength we should 
consider, it is conceivable that the correlation function 



could receive an infinite contribution from them, possi- 
bly diverging. Unless these waves are weighted in inverse 
proportion to their wavelengths, in which case the con- 
tributions to C(9) from larger waves could saturate at 
some finite number. 

Let us analyze this issue closer by invoking the defini- 
tion of the real-space temperature correlation function: 

C(60HAT(To,fLi),AT(ro,n2)) • 

Here, To is our present cosmological time, n^ represents 
the direction of each incoming photon, and 9 is the angle 
between them. For large angles the Sachs- Wolfe effect, 
AT = <I>/3 for adiabatic initial conditions, is the domi- 
nant contribution to C(0). Going to Fourier space and 
using 

($(k)$(q)) = 2ir 2 k- 3 V(k)6 {3) {k + q) , 

as expected from the symmetries of a Gaussian, homo- 
geneous and isotropic universe [6], we find after a bit of 
algebra that 



C{9) 



dk sin kr 



k kr 



V(k), 



(2) 



where r — 2At\/1 — cos 6 and At is the radius of the last 
scattering surface. Essentially all models of the early uni- 
verse predict a primordial power spectrum V{k) which is 
a power law in k. Without loss of generality this function 
can be written as 



V(k) = Ak n ° 



(3) 



where n s is known as the scalar spectral index and A is 
an amplitude to be fixed. Using this power spectrum in 
Eq. (2) we find that after integration that 



C(0) ex (1-cosl 



1 < n, < 3. 



(4) 



There are two important features to be noted in this ex- 
pression. First, the amplitude is a finite function of n s . 
Second, the above correlation function has no root in the 
range 1 < n s < 3 and < 9 < 180° - see Fig. (2). Indeed, 
this agrees with our discussion, since for n s > 1 the power 
spectrum in Eq. (3) filters out low values of k, ensuring 
that arbitrarily large waves will cease to contribute to 
C(9). For a Harrison-Zel'dovich spectrum (n s — 1) the 
above integral has a finite dependence on 9: 



C{9) ex In 



1 



1 - cos 9 



1/2 



(5) 



However, the amplitude here is a logarithmic divergent 
function of k, as one can check by fixing the lower limit 
of the integral in Eq. (2) at fco and carrying a Taylor 
expansion around zero. Evidently, this divergent ampli- 
tude corresponds to a monopole Co, and can be safely 
regularized. The important point here - and also for the 
case n s < 1, although with Eq. (5) replaced by a more 
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Figure 2: Large-angle temperature correlation function for 
different values of the scalar spectral index n s . The curves 
represent analytical expressions of C{9) for a Harrison- 
Zel'dovich (n s = 1, continuous line) and blue (n s = 1.2, 
dashed line) power spectra. The envelopes represent la cos- 
mic variance error bars. These curves were arbitrarily nor- 
malized at 6 — 1° to match the ACDM correlation function 
(thin, dot-dashed line). 



complicated expression - is that the correlation function 
has one root at 9 = 90° (Fig. (2)). Both the infinite 
amplitude and the root of this function agrees with the 
interpretation of V(k) acting as a constant filter which 
allows arbitrarily large waves to contribute to C(9) if 
n s = 1. 

Faced with this filtering property, another important 
question is related to the use we can make of it. Unfortu- 
nately, using the real-space behavior of C{9) to constrain 
n s would be a real challenge. Not only the errors in the 
measurement of C(9) at large angles are strongly cor- 
related, but the observed C{9) is known to suffer from 
a coherent lack of power at these scales [' )] (see also 
Ref. [10] for a recent review of this and other related 
large-angle anomalies.) We are not going to pursue these 
issues here. Instead, we note that the above example 
poses a different but related question: how far beyond 
the radius of the observable universe can a wave affect 
CMB, for a given n s l Since for a blue spectrum (n s > 1) 
arbitrarily large waves will contribute less and less to the 
temperature fluctuation in contrast to a red one (n s < 1), 
we can expect the former contribution to saturate at a 
larger maximum wavelength than the later. Moreover, 
by restricting the amplitude of these waves to the lin- 
ear sector of perturbation theory, this effect provides a 
natural mechanism in terms of which we can put lower 
bounds on the extension of the homogeneous patch of the 
universe. 

The idea of using large (super-Hubble) waves to es- 
timate the radius of the homogeneous universe is not 
new and was first proposed by Grishchuk and Zel'dovich 
(henceforth refereed to as the GZ effect) [11]. In Ref. [12], 
the effect of super-Hubble linear perturbations from a 
scalar field was investigated, where it was found that the 
linear universe is at least 45 times larger than the ob- 



servable universe. In Ref. [13] a single sinusoidal was 
employed, leading to a lower bound on the radius of the 
homogeneous universe of 65 times the present particle 
horizon. Further extensions were pursued in order to 
apply the GZ effect to spatially curved universes [14], 
and in the estimation of possible inflationary remnants in 
CMB [ ]. However, in all of these approaches, and also 
in the original GZ paper, the effect of the scalar spec- 
tral index was completely ignored. Instead, the "wall" 
of the homogeneous universe was assumed to be well 
characterized by a sharp feature in the power spectrum: 
P(k) = AS(k — kcz)- Since the original GZ effect was 
proposed before the formulation of modern inflationary 
models, it is reasonable that one should model a break in 
the scale of homogeneity by a sharp peak in the spectrum. 
From this point of view, the approach we will follow here 
is an natural extension of this idea. Given that inflation 
predicts the existence of filtered stochastic waves, as we 
illustrated above, this very property will dictate the lim- 
its on the size of the homogeneous universe that current 
CMB data can place. In a sense, our proposal is to es- 
timate the radius of the universe by measuring it from 
inside-out, rather than at the edge. 

This work will be organized as follows. In Sec. II we 
will define precisely what we mean by filtered stochas- 
tic waves, and show how they should be normalized on 
arguments of global scale invariance. In Sec. Ill we re- 
strict attention to super-Hubble waves and calculate the 
impact of these waves on the low multipole moments of 
CMB temperature map. Finally, we use the recently re- 
leased Planck data to put constraints on the size of the 
homogeneous universe in Sec. IV. We conclude in Sec. V 
with our final remarks and some further perspectives of 
development. Throughout the text, we use natural units 
in which c = 1 . Our convention for the Fourier transform 
is to democratically distribute factors of (27r) 3 / 2 among 
integrals. 



Note that the dependence on the direction of k appears 
even in an isotropic background setup, although here it 
comes exclusively from the initial conditions ip, and not 
from the spacetime dynamics as would happen, for ex- 
ample, in anisotropic universes [1G, 17]. The function (f> 
can develop a /c-dependence only through the Laplacian 
appearing in the evolution equation [18]: 



<t>" + 3ft (1 + c 2 s )4>' + [2ft' + ft 2 (1 + 3c 2 )] 4> + k 2 c 2 s (j) = 

(6) 
where c 2 is the sound speed of whatever fluids constitute 
the universe, ft is the conformal Hubble parameter and 
the primes mean derivatives with respect to conformal 
time r. The general idea behind inflationary models is 
that inflation has lasted for a finite number of e-folds, 
stretching an initially homogeneous patch of the universe 
of typical size ft -1 up to some radius L. During this 
evolution, the dynamics of <f> can be linked to that of the 
curvature perturbations C,, which is known to be constant 
at super-Hubble scales. Since in this work we shall be 
mainly interested in scales where k <C ft, we will make 
the mildly technical assumption that (j) does not depend 
on k at all, that is, we assume that 



*(T,k) = 0(rMk) 



(7) 



is a good description for the gravitational potential at all 
wavelengths. This is certainly not valid at small scales 
during, let's say, the radiation dominated phase, where 
<j) is an oscillating and decaying function of k. How- 
ever, during inflation, every factor of k in the solution 
of Eq. (6) is accompanied by a factor of 1/a. So, al- 
though we will integrate Eq. (7) over all values of k, it is 
expected that the corrections to the /c-dependence of (p 
are suppressed by the quasi-exponential expansion of the 
inflationary universe [ ]. In conclusion, Eq. (7) ensures 
that $(t, x) is also a separable function of r and x, which 
we can always decompose as 



II. FILTERED STOCHASTIC WAVES 



$(r,x) = </)(t) } w em (x)Y lm (h) . 






(8) 



We will start by considering the standard scenario of 
scalar and adiabatic perturbations evolving in a spatially 
flat, Friedmann-Lemaitre-Roberson- Walker (FLRW) uni- 
verse dominated by cold dark matter, baryons, radia- 
tion and a cosmological constant. In the absence of 
anisotropic stress, the spacetime metric is 

ds 2 = a 2 (r) [-(1 + 2$)dr 2 + (1 - 2^)5 i;j dx i dx j ] 

where $(r, x) is the gauge-invariant gravitational poten- 
tial. Usually, linear perturbation theory is carried in 
Fourier space, where $(t, k) is the relevant quantity. In 
the standard FLRW case, we can split the functional de- 
pendence of this function as a dependence on the initial 
conditions due to inflation, ^(k), times the dependence 
on the dynamical evolution due to gravity, ^(t, k): 

*(T,k) = #T,*)lKk). 



Here, the multipolar coefficients W£ m (x) are understood 
as complex Gaussian random variables satisfying 

(w£ m (x)) = 0, (wtm(x)w£i m i(x)) = |w*(a;)| 5w6 mm i , 

where, it should be noted, the averages are taken at a 
fixed radius x. This ensures that the gravitational poten- 
tial in Eq. (8) is also a Gaussian random variable whose 
two-point correlation function is given by an expression 
similar to Eq. (1). In what follows, we will rewrite the 
multipoles W£ rn (x) as 

^i m (x) = wi(x)(ftm (9) 

where tpi m is a (complex) Gaussian variable satisfying 

{ftmfl'm') = Su'Smm' . (10) 



This decomposition is convenient because it separates 
the gravitational power spectrum wg{x). which is a fixed 
physical quantity, from the purely random phases (pg m . 
This splitting is always possible in a FLRW background 
universe, where CMB is statistically isotropic. 

Our task now is to relate the real-space spectrum wg(x) 
to the scalar spectral index n s . In order to proceed, note 
that if we decompose Eq. (7) similarly to what was done 
in Eqs. (8) and (9), we can relate wi{x) to its Fourier 
dual, zug(k), through the so-called Hankel transform [4]: 



wi{x) 



dk k jg(kx)vu£(k) . 



(11) 



The Fourier-space spectrum zui(k) is related to the ini- 
tial conditions ip(k), which are also Gaussian random 
variables satisfying 

(V>(k)V>(q)) = P(k)S {3) (k + q), P(k) = 2ir 2 Ak n °- i . 

(12) 
Writing a simple expression for V'(k) is virtually impossi- 
ble since, for each wavenumber k, inflation has produced 
different initial conditions with independent phases and 
amplitudes. This is confirmed by the simple observa- 
tion that galaxies tend to cluster differently at different 
scales and in different directions. Since we want to place 
bounds on the size of the homogeneous patch of the uni- 
verse, we can speculate that a typical realization of the 
initial conditions at super-Hubble scales is 



^k = y/Plkje 



ik-L 



(13) 



where L = |L| is the radius of the homogeneous patch. 
Physically, this ansatz represents a standing wave whose 
amplitude is filtered by the square root of the variance 
P(k) - we call it a filtered stochastic wave. We stress that 
Eq. (13) is an ansatz for one realization of ip at a fixed 
wavenumber k, not the full Gaussian field. Nonetheless, 
this ansatz is quite general because we can always decom- 
pose ^(k) as \/P{k) V'(k), and write V'(k) as a Fourier 
series of plane waves e _lk ' L . Evidently, the fact that dif- 
ferent waves have different phases and amplitudes has to 
be taken into account. This can be easily implemented 
in harmonic (rather than in Fourier) space by defining 
wi(k) as the angular coefficients of ipu 1 



y/P(kj< 



-ik-L _ 



^(2^+l)n^(/c)P £ (cos60 



from where it follows that 

w e (k) = {-if^P[k) 3i {kL) 



(14) 



Then, since we are working under the hypothesis of sta- 
tistical isotropy, the randomness of phases and ampli- 
tudes can be implemented by multiplying zui(k) by the 



1 Note that these coefficients can also be interpreted as the angular 
average of Eq. (13): ro^(fc) = n /— i dp ipkPg (m)- 



standard Gaussian variable <pi m . This is achieved auto- 
matically by Eqs. (9) and (11). 

We can now compute the integral in Eq. (11) using 
Eq. (14). Since by hypothesis x <C L, the result can be 
expanded in powers of x/L. We find 



wt{x) = g e (n s )e e + O (e e+2 ) 
where e = x/L <C 1 and 



<)( 



(n.) = 2^-V/ 2 



2+n„ 



r(< + §)r(fcp) 



(15) 



(16) 



This expression, together with Eq. (9) and the solution of 
Eq. (G), completes the construction of the gravitational 
potential at large scales: 

$(t, x) = c/)(t) Y^ 9i(n s )(pi m Yim(ii.) e e . (17) 

As a last comment, note that, in order to obtain Eq. (15), 
we have fixed the amplitude in Eq. (12) as A = L ris+2 . 
This choice ensures that the gravitational potential is a 
homogeneous function of x and L, that is, $(x;L) = 
<f>(\x] XL) for some real A. This is necessary because x 
and L are comoving coordinates, so that their scalings are 
meaningless. Incidentally, this prescription does not fix 
A uniquely, since the choice A = BL Us+2 would also work 
for any constant B. However, this new constant can be 
absorbed in the initial value of 4>(t) entering in Eq. (17), 
which will be properly normalized at CMB decoupling 
later on. 



III. LARGE SCALE TEMPERATURE 
FLUCTUATIONS 

Moving forward, we will now compute the contribution 
of the large scale gravitational potential in Eq. (17) to 
the low multipoles of CMB. The main contribution to 
CMB at large scales can be accounted for through the 
well-known Sachs- Wolfe formula [ ]: 



AT(r ,n) 



1 



$ 



(r,x dcc ) + 



3$(t,x(t)) 
dr 



dr + n • Vul 



(18) 



where 



X ( T ) = (r -r)n. 



Since at CMB decoupling the gravitational potential is 
^(TdocXdcc), we will present our results in powers of 



Qlc 






«i, 



where xa cc = |x(rd oc )|. In order to calculate the mul- 
tipolar decomposition of Eq. (18) we will need to solve 



Eq. (6) numerically, together with the background Fried- 
mann equations. We have chosen a fiducial ACDM model 
with best fit parameters given by Planck data alone [1], 
which are: fl m — 0.3175, z cq — 3402 and z dcc = 1090.43. 
To solve Eq. (6) we use adiabatic initial conditions such 
that 



</>(0) = <f>P, 0'(O)=O, 



(19) 



where p is the primordial value of <j>{r). As we are going 
to see, all contributions to AT/0(Tdec) depend on ratios 
like <J)(t) / <j>(jdcc) which, at large scales, are independent 
of the choice of the value of <fi p . On the other hand, the 
absolute value of cj>(Tdcc) does influence our estimates of 
AT, and so it will have to be fixed. We fix this value by 
restricting our analysis to the linear domain of perturba- 
tion theory at CMB decoupling, that is, we impose the 
condition 



\Sp/pl 



<1. 



(20) 



At large scales, the relation between the energy density 
perturbation and the gravitational potential is 

4TrGa 2 5p sa -3H(<f>' + H0) , 

where we have used the separability condition in Eq. (7) 
to isolate the spatial dependence of Sp. Changing vari- 
ables from r to a and using 3"H 2 = 8TrGpa 2 , the time 
dependence of the fractional energy density becomes 

Sp „ d , ,. 

-^ = -2— M . 

p da 



Defining <j>(a) = </> p f(a) such that /(0) = 1 and /'(0) 

0, we find numerically that \Sp/p\ T 

condition in Eq. (20) then fixes (j){ T dcc) uniquely as 



1.828 <p p . The 



<t>(r dec )= 0.512, (21) 

which, in turn, fixes the ratio L/x dcc , &s we shall see. 



Using Eq. (17), the Sachs- Wolfe temperature fluctuation 
is then written as 



AT sw (t , n) = <j>{T dcc ) ^2 PirnYtmta) S e e dcc , (22) 
where 



(i.ni 



S e = g e (n s )J\f . 



B. Integrated Sachs- Wolfe effect 



(23) 



In order to compute the contribution from the inte- 
grated Sachs- Wolfe (ISW) effect we use the identity 

<& (t, x (t)) = — n • V$ 

OT 

where dot means total derivative with respect to r and 
where we have used dx/dr = — n. Plugging this in the 
ISW expression we readily find 



ATisw(T ,n) = 2$| 



n • V$dr . 



If we align L with the z-axis and write x in spherical 
coordinates, then fi • V = d/dx and we have 

n- V$ = (j){T)^2gnp im Y em (h)— e £ . 



Cm 



Joining all the terms and dropping the monopole 
2$(to,0), we arrive at 



ATisw(r ,n) 
where 



(Tdoc) ^2 PlrnYtmia) %i £dcc (24) 

i.m 



If 



-2g e (n s 



1 



Hr) d 
(r de c) dx \x d , 



dr 



(25) 



C. Doppler effect 



A. Sachs- Wolfe effect 

The Sachs- Wolfe (SW) effect is the simplest contri- 
bution to compute. First, note that at super-Hubble 
scales the radiation plasma has a negligible peculiar 
speed, which implies that S' r — 4$' = [18]. Restrict- 
ing our analysis to adiabatic initial conditions, for which 
<5 r (0,x) = -2$(0,x), it then follows that 



1 



^(TdccXdcc) + $(T dcc ,Xdcc) = $(T dc c,X dec )A/', 



where J\f = 2 — 30(O)/2</>(rdec) needs to be evaluated 
numerically. We find 

Af= ^+0.0644^0.4. 



At the scales we are interested, the mean free path 
of photons is very small, and we can think of baryons 
and radiation as tightly coupled into a single fluid (tight- 
coupling limit). The Doppler effect at direction n then 
involves the baryon-radiation plasma peculiar velocity 

v = V« through fi • Vu|I° . The relation between v 

. . Tdoc . 

and the gravitational potential follows from Einstein's 

equation 



d_ 



(o$) = -4?rGa 3 p(l + u)v 



where p = p m +p r , and u) is the effective equation of state 
of the plasma. Since in our approximation $(r, x) = 
</>(r)?/>(x), the velocity can be rewritten as 

V = /( T ) X! MVlmYtmifye* . 

l,m 



with /(r) = —gp {a<t>) /4irGa 3 p(l + w). Using n • Vv = 
dv/dx, we find 

AT Dop (r , x dcc ) = 0(r dec ) ^ cp em Y im (n)'D i e^ ec , (26) 



where we have introduced 
_ ge(n s ) 



V,= 



Xdec4>{ T dec) 



[f(r )S ei - f(r dec )£} . (27) 



IV. OBSERVATIONAL CONSTRAINTS 

We will now use the recently released Planck data 
to put constraints on the ratio L/xdcc- With the rele- 
vant effects computed, we can gather all contributions to 
AT(to, n) in the simpler expression 



AT(t , n) = 2^ ai m (T )Ye m (n) . 

l,m 

with the multipolar coefficients given by 



i( T o) = 0(T"dc 



(S t 



V, 



(28) 



(29) 



This expression, which is valid up to terms of order 
e dcc ' snows that the stronger constraint on the ratio 
L/x^ec will come from the lower CMB multipoles. Since 
the monopole Co is unobservable, the first contribution 
comes from the intrinsic CMB dipole. However, in what 
concerns super-Hubble waves, the intrinsic dipole pro- 
duced by Eqs. (22) and (24) is exactly canceled, at lower 
order in edec, by the Doppler dipole in Eq. (26). This 
can be seen in the first panel of Fig. (3), where we show 
some plots of the first multipoles ATi entering the total 
temperature fluctuation, and defined by 



A7>(0; n.) = (Se+le + V e )Y eo (d) . 



(30) 



The dipole cancellation is a peculiar feature of the GZ 
effect which happens quite generally in a ACDM universe 
- see Ref. [13] for an explicit demonstration. Here we find 
that it also holds, regardless of n s . As a general behavior, 
we can see that the amplitude of the other multipoles 
grow with growing values of n s . 

It's also instructive to see the individual contributions 
from the SW, ISW and Doppler effects to the multipoles 
in Eq. (30). These are shown separately in each panel of 
Fig. (4). Although not explicit, one can infer from the 
first panel in this figure that the sum of the Doppler and 
ISW contributions exactly cancels the SW term. 

In order to proceed, note that, since ipg m are standard 
Gaussian variables, it follows from Eq. (29) that 

C t = <Xr dc c) 2 fit +%£ + Vef ef ec . 

Using Planck's best fit value n s — 0.9624 and the linear- 
ity prescription of Eq. (21), the quadrupole amplitude 
evaluates numerically to 



C 2 = 0.067 



\17) 



(31) 
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Figure 3: Total multipolar contribution from super-Hubble 
waves to the CMB as a function of 6. From top to bottom 
we show the dipole ATi, quadrupole AT2 and octopole AT3. 
Continuous, dashed and dotted curves represent the values 
n s = (0.5, 1.0, 1.5), respectively. The dipole is exactly zero in 
a ACDM universe; see the text for details. 



The nominal value of the best fit quadrupole measured 
by Planck is 



2(2 + 1) 
2tt 



C 2 = 299.495 x 10~ 12 . 



In the (unlikely) event that the full quadrupole is at- 
tributable to a super-Hubble wave as defined here, we 
obtain as a lower bound L/xdcc ^ 120. Any fraction of 
the measured quadrupole smaller than one would only 
raise this value, so this is actually a conservative lower 
bound. In fact, since this ratio decreases with increas- 
ing C2, a more conservative bound would be given by 
the measured C'2 with error bars. The Planck team has 
shown that, at 68% confidence level, the quadrupole is 



2(2 + 1] 

2tt 



C 2 = (299.495 tlUi^)x 10" 



Taking the upper error bar, we find 



L 



>87. 



(32) 



Just for comparison, let us see what would a purely theo- 
retical estimation look like. Using CAMB [21] to evaluate 
the quadrupole for the fiducial parameters used in this 
work, we find 2(2 + 1)C 2 /2tt ~ 1136 x 10" 12 . This gives 
L/x^cc ^ 86, which is essentially the same as Eq. (32). 
Adding la cosmic variance error bar to this value will 
lower it to L/xdcc ^ 76, which is again not too far from 
the observational bound. Interestingly, this shows that 
the tension between the measured and predicted value 
of the quadrupole Ci is not so drastic when viewed from 








Figure 4: Individual multipolar contribution from the SW 
(continuous), ISW (dotted) and Doppler (dashed) effects as 
functions of 9. A simple inspection shows that the sum of 
the curves in the first panel gives zero. See the first panel 
in Fig. (3). These plots were generated for a Harrison- 
Zel'dovich spectrum and show, from top to bottom, the 
dipole, quadrupole and octopole of each term in Eq. (30). 




Figure 5: Quadrupolar bound on the ratio L/xdcc as a func- 
tion of n s based on Planck data with lcr (lower) error bar (see 
the text for details). Models with red spectrum (n s < 1) are 
more sensible to super-Hubble waves, implying in a smaller 
lower limit to the size of the universe. Blue spectrum (n s > 1) 
models give less weight to larger waves, implying in a higher 
lower bound. Vertical lines represent the Harrison-Zel'dovich 
(continuous) and Planck's best fit (dashed) spectra. The 
lower dashed line represent the ACDM model with lcr (lower) 
cosmic variance error bar. 



FINAL REMARKS 



the real space perspective. This sort of analysis reinforces 
the idea that a real space approach to CMB [3, 4, 10] de- 
serves attention, even if in principle it contains the same 
information as found in the more popular harmonic ap- 
proach. 

What does the bound in Eq. (32) have to say regard- 
ing the duration of inflation? Supposing that inflation 
has lasted N m [ n = lnaf/a; e-foldings, just enough to solve 
the horizon and flatness problems, the size of the homoge- 
neous patch larger than the Hubble horizon at the onset 
of inflation would be only slightly larger than our hori- 
zon today. Therefore, any departure from the lower limit 
L/xdcc ^ 1 would imply in a lower bound to the num- 
ber of e-foldings, that is N > N min + In 87 w N min + 4. 
Incidentally, this suggests that CMB data is barely at 
the observational window for signatures of primordial 
anisotropics [22]. 

Finally, going back to the original question raised in 
this work, we can check that the bounds on the size of 
the homogeneous universe is a function of the scalar spec- 
tral index. Since red spectra are more sensitive to larger 
waves, the CMB quadrupole saturates at smaller super- 
Hubble lengths than the saturation length of the scale in- 
variant (Harrison-Zel'dovich) spectrum. Blue spectra, on 
the other hand, are less sensitive to larger waves, making 
the measured quadrupole value saturate at larger super- 
Hubble lengths. This behavior is depicted in Fig. (5). 



The primordial power spectrum predicted by most in- 
flationary models displays a filtering property noticeable 
at very large cosmological scales. This can be demon- 
strated by writing the explicit dependence of the angular 
two-point correlation function with the scalar spectral 
index n s . By employing the Sachs- Wolfe effect, we have 
shown that this filtering property produces a correlation 
function which is always positive if n s > 1, whereas large 
scale correlations will unavoidably change sign if n s < 1. 
This mechanism implies that the spectral index plays 
an important role in estimating the size of the homo- 
geneous patch of the universe through, for example, the 
Grishchuk-Zel'dovich effect, where the size of the uni- 
verse is estimated through the modulation it introduces 
at low CMB multipoles. In this work we have extended 
the GZ effect to comply with stochastic, super-Hubble 
waves. Because these waves have amplitudes which are 
filtered by the primordial power spectrum, the estimates 
on the allowed size of the largest super-Hubble wave de- 
pend on the value of the scalar spectral index n s . Infla- 
tionary models predicting an index n s < 1 have a spec- 
trum which is more sensible to super-Hubble waves as 
compared to models predicting n s > 1. This means that 
the contribution of red waves to the low-£ CMB multi- 
poles will saturate at shorter super-Hubble length than 
the contribution of blue ones, implying in lower bounds to 
the size of the homogeneous patch of the universe. Using 
Planck's best fit value for n s and the quadrupole C2, we 
have found that the homogeneous universe extends to a 
distance at least 87 times bigger than our present optical 
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horizon. This result is compatible with analytical esti- 
mates found in previous works [12, 13, 15], and one order 
of magnitude smaller than the statistical analysis car- 
ried in [23], where the standard GZ effect was employed. 
However, we emphasize that our approach is based on a 
different philosophy. Instead of evoking a phenomenolog- 
ical sharp feature in the primordial power spectrum, we 
use the very inflationary power spectrum as cosmologi- 
cally motivated spectral filter. In other words, we have 
measured the size of the universe from its interior up to a 
saturation quadrupole value fixed by CMB data. As one 
last comment, we would like to point to some possible 
applications of the results presented here. Recently, the 
authors of Ref. [24] claimed that frequency distortions 
on the blackbody spectrum of CMB, such as those pro- 
duced by the kinetic Sunyaev-Zel'dovich effect and the 
Compton-j/ distortion, could be used as complementary 



data to the GZ effect. Indeed, in the light of our results, 
complementary information on the low-£ properties of the 
universe could used to independently constrain the scalar 
spectral index, thus narrowing the space of cosmological 
parameters yet in the regime of Gaussian and isotropic 
perturbations. 
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